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ABSTRACT 


With the increased rate of urbanization and industrialization, the Nigerian 
electric power system is being put under pressure, high power losses which 
has led to fluctuation in voltage level. In this work the Nigeria 330Kv network 
of 42 bus system was considered. Data for the analysis were obtained from 
Transmission Company of Nigeria / National Control Centre, Osogbo (TCN) 
and MATLAB/PSAT software with newton Raphson's solution method 
embedded in it was used to carry out the analysis. The results of the analysis 
showed that many of the bus voltages were outside the voltage limits of ±5% 
i.e 0.95pu- 1.05pu. After compensation with Flexible AC Transmission System 
devices (FACTs) namely Unified Power Flow Controller (UPFC) which can be 
used to control power flow on a transmission line, the voltage profile almost 
flat with bus voltages within acceptable voltage limits. It is also evident that 
the UPFC is device that can used to combat the voltage problem in the Nigerian 
electric power system. 



© 


1. INTRODUCTION 

Electric power generation and transmission is a complex process, wherever 
power is to be transferred, the two main components are active and reactive 
power. In a three phase ac power system active and reactive power flows from 
the generating station to the load centers through different transmission lines 
and networks buses. The active and reactive power flow in transmission line is 
called power flow or load flow. The Nigerian power system is characterized by 
high power loss that affects the reliability and quality of power supply. 
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Today urbanization and rapid industrialization are putting 
power utilities under pressure. The consumers at every level 
including domestic, commercial and industrial users have 
been on the receiving end of the voltage instability problem 
in Nigeria electricity grid, experiencing major discomforts in 
operating their electrical and electronic devices as well as 
unable to run smoothly their day-to-day activities since they 
depend largely on stable and quality electricity supply. 

Despite the reform exercise and other corrective measures 
put in place by the Federal Government of Nigeria to address 
the operational challenges facing her electricity grid, voltage 
instability problem still persists. Most of voltage instability 
problems in the Nigeria electricity grid can be traced to 
inadequate and inefficient reactive power compensation in 
the system. Therefore, the need to devise means of 
compensating reactive power in Nigeria electricity supply 
network becomes highly imperative for efficient operation of 
the system. To address this, adoption of advanced control 
technologies such as FACTS is an important option for the 
Nigeria electricity system reactive power compensation in 
view of the slow response of the conventional traditional 
methods for improving power system performance [1,2]. 

FACTS technology has a lots of benefits, such as greater 
power flow control ability, increased in the loading of 
existing transmission circuits, damping of power system 
oscillations, has the less cost than other alternative 
techniques of transmission system is used [3,4,5,6,7,8]. The 


ability of FACTS to control the line impedance and the nodal 
voltage magnitudes and phase angles at both the sending 
and the receiving ends of key transmission lines, with almost 
no delay, has significantly increased the transmission 
capabilities of the network while considerably enhancing the 
security of the system. Some of the common FACTS 
controllers in use are Interphase Power Controller (IPC), 
Static Synchronous Compensator (STATCOM), Static Var 
Compensator (SVC), Thyristor Controlled Breaking Reactor 
(TCR), Thyristor Controlled Series Capacitor (TCSC), Static 
Synchronous Series Compensator (SSSC), Thyristor Switched 
Series Reactor (TSSR), Unified Power Flow Controller 
(UPFC). These controllers have their distinct characteristics 
which are suitable for different applications and due to these 
various researches have been carried out to examine their 
potential applications and benefits [9,10,11,12,13,14,15, 
16], 

Lots of works have been done to improve the voltage profile 
of the Nigerian electric power using other conventional 
methods of compensation, FACTS devices like STATCOM, 
SSSC, SVC but in this study, our aim is to apply UPFC for 
improvement of voltage profile in the Nigerian power 
system. 

2. Unified Power Flow Control (UPFC) 

Unified Power Flow Controller (UPFC) is power electronics 
based system that can provide the control of the 
transmission line impedance, phase angle and reactive 
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power. This versatility of the UPFC makes it a prime FACTS 
device that can provide many of the control functions 
required to solve a wide range of dynamic and steady state 
problems encountered in power systems. Combining the 
STATCOM and the SSSC into a single device with a common 
control system represents the third generation of FACTS 
known as Unified Power Flow Controller (UPFC). It has the 
unique ability to control real and reactive power flow 
independently. The basic operation principle diagram of the 
UPFC is shown in Figure 1, and has been described in open 
literature [17, 18, 19, 20, 21, 22]. The Unified Power Flow 
Controller (UPFC) made out of two Voltage Source 
Converters (VSCs), which are connected through a common 
DC link capacitor. 


From figure 1, the basic function of converter 1 is to supply 
or absorb the real power demanded by converter 2 at the 
common dc link. This dc link power is converted back to ac 
and coupled to the transmission line via a shunt-connected 
transformer. Converter 1 can also generate or absorb 
controllable reactive power, if it is desired, and there by it 
can provide independent shunt reactive compensation for 
the line. It is important to note that whereas there is a closed 
"direct" path for the real power negotiated by the action of 
series voltage injection through Inverters 1 and 2 back to the 
line, the corresponding reactive power exchanged is 
supplied or absorbed locally by converter 2 and therefore it 
does not flow through the line. Thus, converter 1 can be 
operated at a unity power factor or be controlled to have a 
reactive power exchange with the line independently of the 
reactive power exchanged by the by the converter 2. This 
means there is no continuous reactive power flow through 
UPFC. 


Shunt Side 


Series 

transformer 


Series Side 



Figurel: Implementation of the UPFC by back to back 
source converters. 


2.1. Mode Of Operation 


Operation of the UPFC from the standpoint of conventional 
power transmission based on reactive shunt compensation, 
series compensation, and phase shifting, the UPFC can fulfill 
these functions and thereby meet multiple control objectives 
by adding the injected voltage V pq , with appropriate 
amplitude and phase angle, to the terminal voltage V. Using 
phasor representation, the basic UPFC power flow control 
functions are illustrated in Figure 2 



Figure2: Phasor diagram representation of modes of 
operation of UPFC 

A. Terminal Voltage Regulation Terminal voltage 
regulation, similar to that obtainable with a transformer 
tap- changer having infinitely small steps. Where 
V pq = AV is injected in phase or antiphase with V. 

B. Line impedance compensation or Series capacitor 
compensation V pq = V c is in quadrature with the line 
current. 

C. Transmission angle regulation, (phase shifting) V pq = V a 
is injected with angular relationship with respect to V 
that achieves the desired phase shift (advance or retard) 
without any change in magnitude. 

D. Simultaneous control of voltage, impedance, and angle 
Multifunction power flow control, executed by 
simultaneous terminal voltage regulation, series 
capacitive compensation, and phase shifting where 
V pq = AV + V C + V a . 

3. Problem Formulation 

Due to the complexity that arises in solving the equations, 
since the functions for real and reactive powers are 
expressed in terms of non-linear algebraic equations. 
Iterative methods for solving equations shall be used in 
solving load flow problems as related to this work. A 
Newton-Raphson iterative technique was of load flow was 
used to simulate and investigate the power and voltage at 
each bus. A MATLAB/PSAT software was also used to run 
the simulation as it has the Newton-Raphson technique 
embedded in it. 

3.1 Modelling of UPFC for Active and Reactive Power 
Evaluation 


The modelling is carried out for the evaluation of the real 
and reactive power using the equivalent circuit of figure 3 



Figure3. Equivalent circuit of a UPFC between two 
busesiand j 


The active and reactive powers of the series convertor 
(SSSC) are as follow. 

Sse = Pse + jQse = V se l*j=V se [Y* n V\ + VjtVj 
+ Yij7* se ].(1) 

Expanding equation 1 and separating the real and imaginary 
points we have that 

Pse = /V S e/ 2 Gij + /v se / /Vi [Gy cos( S se - 5;) + BjiSin (S se 
- Sse)] +/V se //V j /[G ij (5 se - 5j) 

+ Bjj sin(5 se - Sj) ] ... (2) 


Qse=-/Vse/ 2 /B i j/+/V se //V i /[G ji sin(<? se - Si) - B n cos(5 se - 
Si+/Vse//VJ/[GJis i n Sse- Sj+BjicosSse- Sj] 

.( 3 ) 
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Also, the active and reactive power of the shunt controller 
(STATCOM] are obtained as 

Ssh = Psh + jQsh = V sh r sh = -v sh v; h [v; h Y s \] .[4] 

Expanding equation 4 and separating the real and imaginary 
parts we have 
Psh — 

~/Vsfi/2GiO+ /Vsh//V//6’/0cos(8sh-8i+B/0s/n (5sh-Sj)J 

.(5) 

Q sll =-/V sh / 2 /B i0 /+/V sh //V i /[G i0 sin (S sh - S t ) - 

Bi o cos (8 sh - 5j)] .... (6) 

Since we assume lossless converters, the UPFC neither 
aborts nor injects active power with respect to the AC 
system. Hence the constraint equation is 
Psh + Pse = 0 . C 7 ) 

Where 

P se is the series converter real power flow 

P sh is the Shunt converter real power flow 

Q se is the Series converter reactive power flow 

Q sh is the Shunt converter reactive power flow 

V se is the Injected series voltage source 

V sh is the Injected shut voltage source 

I sh is the Current flowing through the UPFC shut converter 

G i0 is the self conductance of the shunt converter 

B io is the self susceptance of the shunt converter 

8 se is the Voltage angle of the injected series voltage source 

8 sh is the voltage angle of the injected shunt voltage source 

3.2 Load Flow Analysis with Newton-Raphson Method 

Load flow studies are one of the most important aspects of 
power system planning and operation. The load flow gives us 
the sinusoidal steady state of the entire system - voltages, 
real and reactive power generated and absorbed and line 
losses. Newton-Raphson's solution method was used to carry 
out the analysis because of its sparsity, fast convergence and 
simplicity attribute as compared to other solution methods. 
Consider the line diagram of a two bus system shown in 
figure 4 


busi 

vi 




JTO_H 

P.i+jQj v i 


Figure4: Two bus system line diagram o transmission 

line 


The power-flow problem discusse in this section will be 
presented in terms of the Y bus Ybus matrix whose elements 
are of the form 

Yij = I YtjW* 11 = I YijIcosQij + j | Yy|sin 0 y = Gy + 

JBij .(8) 


Let Pi and Qi denote the net real andreactive power entering 
the network at bus i. Then the complex conjugate of the 


power injected at bus is given b Si = + jQi.(11] 

Pi- jQi= v 1 *z£ =1 Y in v n = v 1 r L .(12) 

Substituting equation 8 and 9 in 12 e have 

Pi - jQi £n=l/Yin //v n /< e in + 5 n - 5,.(13) 

From the preceding equation we obtain the following form of 
the power-flow equations: 

Pi =Z2=l \Yin ViV n cos (6 in + d n - Si) . (14] 

Qi =Zn =1 I Yin ViVnlsin (0 in + S n - S t ) . (15] 


Where 

Vi is the voltage at bus i 
Vj is the voltage at bus j 
Si is the voltage angle at bus i 
Sj is the voltage angle at bus j 

Expanding equation 13, we have 

Pi -JQi = V t G tJ - [Videos (Si - 8j) 

— VjVj sin (Si — 6j)]G tJ + V?jB tj 
- [Videos (S f - S j) 

+ ViVjj sin (Si - 8j)\jBij 

= V, 2 Gjj - ViVjGijCos (Si - S j) + ViVj[GijSin(8i - 8j) + 

VfjBtj - V i V j jB ij cos(8 i - 8j ) + V^sin^ - 8 } ) .(16) 

Separating the real and imaginary parts from both sending 
and receiving end we have that: 

For sending Bus 

Pij = ViGij - VtVjiGijCos (Si - 8j)+ By sin (6 ( - 5,)] 

.(17) 

Qij = -V?Bij - ViVjlGijSin (5, - 5 ; ) - By cos(5i - 5,)] 

.(18) 


Similarly, for the receiving bus 

Pji = Vj 2 Gji ~ ViVjlGjt cos(8j - 6;) - B Jt cos(S ; - Si)] 

.(19) 

Qjt = -VfBjt ~ Vji[Gji cos (8j - Si) - By cos(Si - S,)] 

.( 20 ) 

The objective of the Newton-Raphson method is to produce 
values for | V,| and Si that will match the prescribed Prf/and Qdi 
as determined from Equations (14] and (15]. At each 
iteration of the method, new estimates of \Vi\ and Si for the 
non-slack buses (/= 2, 3, ..., N ] are generated. At the end of 
each iteration, the power mismatch is given by 


APi = Pi,sch - Pu .(21] 

AQi-Qi,sc\v ■ Qi .(22] 


For i, k = 1, 2,..., N. Let the voltae at bus I be denoted by 

Vi = I V;|e ;5i = |V;|(cosS; + jsind^ .(9] 

For i= 1, 2,..., N 

The net current injected into the etwork at bus I in terms of 

the elements Tin of the Y bus is determined by 

Ii= Y u Y 1 + Y I2 Y 2 + ...+ Y iN Y N Y,k=i Yj n ,i=l, 2...,n.(10) 


The complex power flow equations for uncompensated 
transmission system in solved by Newton -Raphson's 
iterative method are defined thus. 

Pn= Pi - Pdi = Y.U/VJ/ VJ/Y in /Cos (0 ln + 6 n - 
6j) i = 1,2..n.(23) 
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Qn = Qi ~ Pdi = Zn=i/Vi// VJ/Ym/CoS (0 in + 6 n - 
5j) i = 1,2,..n... (24) 

For compensated transmission system, equation 24 is 
modified and is defined by 

Qn = Qi - Qdi + Q * =~ Zn=i/Vi//V n //Y in /sin (0j„ + 

6„ - 6*).(25) 


Y in = element of the bus admittance matrix 
S; = Voltage angle at bus i 
S n = Voltage angle at nth Bus 
0; n = Angle associated with Y in 

Then the new estimates for the bus voltages are 

?(fc) 


Where 

P t = Real power generated at the i th bus 
Qi= Reactive power generated at the i th Bus 
?di = Real power consumed at the i th bus 
V n = Voltage of the n th bus. 

Q di = Reactive power consumed at the i th bus 
Vi = Voltage at bus i 


(26) 


5 , (k_1) = S, (k) + AS 

l l L 

\Vi k+1 \ = \v[ D) \ + A|V i (k) |.(27) 


The solution of equation enables us to identify weak points 
in the system where the voltages magnitude lies outside the 
voltage limits of ±5%. 


3.3. Modelling load Flow Equations with UPFC 

The linearized system of power flow equations for UPFC in connection with the rest of the network in obtained using Newton- 
Raphson power flow equation given by 
[F(x)] = [J][AX] .(28) 


Where 

F (x) is the Power and control mismatch vectors 

Ax is the Incremental vector of state variables Where AX is the solution vector and it is given by 


AX= A Sf 

r 




V. 


J 


[FX] = 


AR 


r 


A6 j 


AP; 

A/V sft 


A Qj 

A /Vj/ 

.29 

AQ; 

A 5 se 


AP;; 

< 


AQ;i 

A5 sft 


AP 


'N 


V 




. 30 


J is the matrix of partial derivative of F(x) with respect to x, 

Let us assume the UPFC is connected to node i and the power and the power system is connected to node j. UPFC is required to 
control voltage magnitude of node i and active power low from node j to node i. Reactive power is infected at node j. Here we 
can re-write equation 28 as 


A Pi 
A Pj 

&Qi = 

A Qj 
A Pji 
*Qjt 
AP 



23_ 

d/V sh / 


dQ i dQ t 
dbj d/V sh / 


dQ j 
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4. SIMULATION AND RESULT ANALYSIS 

Overview of Nigeria 330kV Transmission Network used for the case study 

The input data for the power flow analysis include the bus data that is real and reactive powers of the generator buses, 
transmission line data (impedance of lines), voltages and transformer/load data obtained from Power Holding Company of 
Nigeria (PHCN) are as presented in Tables contained in the appendix. They are used to carry out the analysis. The single-line 
diagram of the existing 330kV Nigeria transmission network used as the case study is as shown in Figure 5. It has 42 buses with 
nine generating station. The Egbin power station was chosen as the slack bus because it has the highest generating capacity of 
1320MW. 

4.1 SIMULATION STEPS 

MATLAB/PSAT software was used to simulate the proposed model. The model of simulation of the proposes model using the 
developed MATLAB code are as follows. The algorithm is represented in figure 4 

> Run load flow on the proposed model before and after incorporation with UPFC. 

> Plot the graph of voltages (p.u) against the buses. 

> Compare the result in graph. 



Figure4: Flow Chart for the Analysis of UPFC Compensation Algorithm. 



Figure5: Modelled Nigerian 330kV power system on PSAT 
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5. RESULTS 

Load flow analysis was carried out on the network to investigate the buses outside voltage limits and also losses along the lines. 
Results of the simulations are in figures 5, 6, 7 and table 1. 


Tablel: Bus voltages and angles with and without UPFC 


Bus 

Number 

Bus Name 

Voltage (p.u) 
without UPFC 

Voltage (p.u) 
with UPFC 

Voltage angle 
without UPFC 

Voltage angle 
with UPFC 

Busl 

B. Kebbi 

0.92327917 

0.98177 

-0.26392 

-0.12408 

BuslO 

Ganmo 

0.993467034 

0.99455 

-0.13674 

-0.12508 

Busll 

Mando 

0.936016761 

0.989 

-0.40837 

-0.37615 

Busl2 

Katampe 

0.96932409 

0.97002 

-0.30952 

-0.28767 

Busl3 

Gwagwalada 

0.972349829 

0.9732 

-0.29775 

-0.27645 

Busl4 

Olorunsogo 

1 

1 

-0.09085 

-0.08538 

Busl5 

Akangba 

0.967384227 

0.98276 

-0.10562 

-0.10388 

Busl6 

Egbin 

1.012 

1.012 

0 

0 

Busl7 

Omotosho 

1 

1 

-0.0174 

-0.01423 

Busl8 

Oke-Aro 

0.980006451 

1 

-0.07579 

-0.07584 

Busl9 

Benin 

1.01126355 

1.0143 

-0.01387 

-0.01021 

Bus2 

Kainji 

1 

1 

-0.10788 

-0.09377 

Bus20 

Kano 

0.912025045 

0.98415 

-0.42084 

-0.3865 

Bus21 

Jos 

0.876609365 

0.98555 

-0.47534 

-0.42692 

Bus22 

Lokoja 

0.967875943 

0.96994 

-0.21473 

-0.1999 

Bus23 

Aja 

1.00915049 

1.0092 

-0.00437 

-0.00437 

Bus24 

Onitsha 

0.989343263 

1.0043 

0.006314 

0.005664 

Bus25 

Ajaokuta 

0.97281235 

0.97499 

-0.18027 

-0.16727 

Bus26 

Delta 

1.012 

1.012 

0.009153 

0.013138 

Bus27 

Sapele 

1.012 

1.012 

0.000476 

0.004477 

Bus28 

Makurdi 

0.893989619 

0.98314 

-0.29127 

-0.26873 

Bus29 

Gombe 

0.812150442 

0.99 

-0.71929 

-0.61011 

Bus3 

Jebba 

0.999264259 

0.99943 

-0.13043 

-0.11731 

Bus30 

New Haven 

0.916384533 

0.98 

-0.10314 

-0.10298 

Bus31 

Okpai 

1.012 

1.012 

0.044018 

0.044556 

Bus32 

Alaoji 

1 

1 

0.172816 

0.168112 

Bus33 

Geregu 

1 

1 

-0.1154 

-0.10237 

Bus34 

Aladji 

1.007415339 

1.0074 

-0.00267 

0.00132 

Bus35 

Ugwuaji 

0.914546079 

0.97883 

-0.10913 

-0.1066 

Bus36 

Yola 

0.801049641 

0.98458 

-0.75828 

-0.63654 

Bus37 

Damaturu 

0.804956113 

0.9833 

-0.75371 

-0.62535 

Bus38 

Afam 

1.003 

1.003 

0.196087 

0.191143 

Bus39 

IkotEkpene 

0.980610751 

0.9889 

0.145729 

0.139584 

Bus4 

Jebba GS 

1 

1 

-0.12617 

-0.11303 

Bus40 

Adiabor 

0.986716932 

0.98844 

0.19726 

0.191375 

Bus41 

Odukpani 

0.991 

0.991 

0.214988 

0.209277 

Bus42 

Sakete 

0.954470925 

0.97019 

-0.13171 

-0.12916 

Bus5 

Shiroro 

1 

1 

-0.27781 

-0.25354 

Bus6 

Osogbo 

0.996650962 

0.99881 

-0.11878 

-0.109 

Bus7 

Aeyede 

0.972034158 

0.97309 

-0.12932 

-0.12185 

Bus8 

Ikeja West 

0.973856683 

0.98913 

-0.09482 

-0.09342 

Bus9 

Ihovbor 

1 

1 

0.049304 

.056249 


1.2 

1 

„ 0.8 

O 0.6 
> 

2 0.4 

0.2 

0 


-m TTOtm - ip rc 
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Figure6: Bar plot of bus voltage with and without UPFC 
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6. DISCUSSION 

The analysis of Nigeria 330kV 42 bus network using 
Newton- Raphson's power flow solution algorithm with 
MATLAB/PSAT software was successfully completed. The 
results obtained revealed the weak buses with values 
outside the statutory limit of 0.95p.u. (313.5kV) and 1.05p.u. 
(346.kV). Four (4] UPFC devices where placed at some weak 
buses in the network namely Birnin Kebbi, Kano, New haven 
and Gombe which resulted to almost a flat voltage profile 
with all the buses within voltage limits as in figure 6. The bar 
plot in the figure 6, shows the comparism of the bus voltages 
before and after compensation by UPFC while figure 7 is the 
corresponding voltage angle of the bus voltages before and 
after compensation. The voltage profile of the system was 
improved as a result of the series compensation voltages 
added to the system with controllable magnitude by the 
UPFC. 
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7. CONCLUSION 

The Nigerian 330kV transmission system associated with 
various challenges like instability of the system as a result of 
voltage profile violation, transmission line inefficiency, 
problem of long transmission lines, network being stretched 
beyond thermal limit, and poor power quality that causes 
constant power failure in Nigeria power system were 
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